Probing MeV Dark Matter at Low— Energy e'^e Colliders 
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It has been suggested that the pair annihilation of Dark Matter particles x with mass between 0.5 
and 20 MeV into e^e~ pairs could be responsible for the excess flux (detected by the INTEGRAL 
satellite) of 511 keV photons coming from the central region of our galaxy. The simplest way to 
achieve the required cross section while respecting existing constraints is to introduce a new vector 
boson U with mass Mu below a few hundred MeV. We point out that over most of the allowed 
parameter space, the process e"^e~ — > Uy, followed by the decay of U into either an e"^e~ pair or an 
invisible {vi' or xx) channel, should lead to signals detectable by current _B— factory experiments. A 
smaller, but still substantial, region of parameter space can also be probed at the $ factory DA"I>NE. 

PACS numbers: 12.60.-i,13.66.Hk,95.35.+d 



Within the context of Einsteinian gravity, evidence 
for the existence of Dark Matter (DM) is overwhelm- 
ing I']. Analyses of the cosmic microwave background 
anisotropy, and other data on the large scale structure 
of the Universe, have determined many cosmological pa- 
rameters with unprecedented precision Along with 
analyses of Big Bang nucleosynthesis (BBN) 3], this data 
also show the DM must be largely non-baryonic. Since 
neutrinos can contribute only a small fraction 0, this 
strongly points towards the existence of an exotic, neu- 
tral, stable particle x, with relic density fl^ satisfying 2] 



0.113 ± 0.0085 (atlcr). 



(1) 



where the scaled Hubble constant ft, ~ 0.7 2]. 

Dark Matter particles should clump due to gravita- 
tional attraction. At the galactic center, their density 
might be so high that their annihilation into lighter, 
known particles could lead to visible signals. Final states 
containing hard photons play a special role in this, since 
photons travel in straight lines and are easy to detect. 

Looking for an excess of hard photons from the cen- 
ter of our galaxy, the INTEGRAL satellite indeed ob- 
served a large flux of photons with energy of 511 keV 
01 . This sharp line can only come from the annihilation 
of non-relativistic e~^e~ pairs. However, most estimates 
of positron production by astrophysical sources fall well 
short of the required flux giving rise to the speculation 
that the annihilation of light DM particles x into e+e" 
final states could be responsible for this signal. While the 
mass TOj^ can exceed substantially (since the positrons 
produced in x annihilation quickly lose energy through 
scattering on neutral atoms), it has been pointed out |h| 
that annihilation into e"'"e~7 final states, in spite of being 
a higher-order process, would over-produce MeV photons 
if > 20 MeV, leading to 



<m^<20 MeV . 



(2) 

Recently, it has been argued that there exists evidence 
for a non-astrophysical source of MeV photons, in which 
case the upper end of the range (j^J would be favored. 



In order to produce the required flux of positrons, one 
needs an annihilation cross section (7ann(xX ~^ e+e^) 



10"-^ fb < WCTa 



(to^/1 MeV)" 



K < 1 fb. 



(3) 



Here v is the relative velocity of the two x particles in 
their cms frame, and a is to be computed at v ^ 10~"^c. 
Note that models of the galaxy fix the DM mass density; 
the number density, which enters quadratically in the 
calculation of the positron flux, therefore scales as 'm~^. 
Finally, k = 1 if x is self-conjugate (i.e. a Majorana 
particle), whereas k = 2 if x 7^ Xj since then only half of 
all encounters of DM particles can lead to annihilation 
events. To be on the safe side, in Q we have expanded 
the range given in the original publication by an order 
of magnitude in either direction. This may be overly 
conservative , since one here only needs the DM density 
averaged over a significant volume, which is thought to 
be better known than that right at the galactic center. 

A cross section in the range ^ implies that e~^e~ ^ 
XX reactions were in equilibrium down to temperatures 
well below m^. Since successful BBN requires a starting 
temperature T > 0.7 MeV 0, it is safe to assume that x 
particles indeed were in thermal equilibrium. Their relic 
density then turns out to be 0| inversely proportional 
to the thermal average of their total annihilation cross 
section into all final states containing only SM particles: 



(4) 



Given the mass constraint to the leading order, an- 
nihilation is possible only into e+e~ and final states; 
the first channel must exist, since x particles are sup- 
posed to annihilate into e"'"e~ pairs even today. Assum- 
ing X forms (nearly) all DM, the relic density resulting 
from eq.Q) must fall in the range Q). This constraint, 
interpreted using eq.Q), and (jSjl are compatible only if 
the present aann is strongly suppressed compared to that 
at decoupling. This is most easily achieved 0, 
if XX annihilation proceeds only from a P— wave initial 
state, in which case v CTann oc v^; note that ~ 0.1 when 
X particles decoupled, while ^ 10^^ today. 



2 



In a renormalizable theory, some particle must mediate 
XX e+e~ annihilation. The simplest possibility 

urn 

is to introduce a light spin-1 boson U coupling to both 
e+e~ and XX states. If x is a Majorana spin-1/2 fermion 
or complex scalar, cr^nnixX ^ //) is given by 
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where — 1 (0) for spinor (scalar), /?/ = y/l — Arn?^ / s, 
n/ = and S/ = gj^ + gj^, with g/^ and gj^ 

being the left- and right-handed Uff couplings. The 
Uxx coupling g^ is purely axial vector for a Majorana x- 
The first line in eq.jS]) is a pure P— wave contribution. We 
discount, henceforth, a Dirac x as it would, in general, 
have a large 5'— wave contribution to a thereby making 
it difficult to reconcile the constraints IpJ and ©. 

Since Tu <C Mu for realistic coupHngs, the usual [13 
non-relativistic expansion of the cross section in powers 
of V breaks down if 2mx — Mjj and the thermal averaging 
has to be done numerically 13]. The strong velocity de- 
pendence of the cross section (s — 4m^ ~ u^m^ if ^ 1) 
in eq.Q also has to be treated properly. For simplicity, 
we assume a thermal velocity distribution, with a rather 



high temperature T — 10 



smaller temperatures, 



which are probably more realistic, would require larger 
couplings, which would be easier to test at e^e" colliders. 
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FIG. 1: Parameter space for a Majorana x with g^j^ — — 
and g^ — lOjjej;. In between the solid curves, x has the 
correct relic density, and the correct cross section to explain 
the flux of 511 keV photons emerging from the galactic center; 
the red (black) curves correspond to 2171^ < {>)Mu- The 
dotted (blue) line indicates the upper bound on (je^ from 
(ge — 2) measurements. The dashed curves show the maximal 
sensitivity of DA<E>NE to e^e~ Uj production, for U 
e+e~ (upper, dark green, curve) and U decaying invisibly 
(lower, magenta curve). The DM constraints are essentially 
independent of the ratio gx/9eRj whereas the ge — 2 constraint 
as well as the sensitivity limits are independent of g^ ■ Results 
for = Sej, are similar to those for scalar x (Fig- 2). 



Figs. 1 (2) show the parameter space spanned by Mu 
and the product of couplings gxdcR for the case that x 
is a Majorana fermion (complex scalar). The regions al- 
lowed by the constraints (|2JEll are enclosed by the solid 
curves, with distinct allowed regions for 2m^ < Mjj 
(2mp(, > Mjj). We consider only ge^ = since this im- 
plies g„^ = if [/ is a gauge boson of a gauge group 
Gu that simply multiplies the gauge group of the Stan- 
dard Model (SM). Note that ve scattering data imposes 
the very strong constraint [l^. g„^ \/ 5ei + 9eR < M^jCp- 
For (7ei = and gej^ = 0, this would exclude the entire 
DM-allowed range (which is invariant under g^j^ ^ gcL ) 
of Fig. 1, and most of Fig. 2. Finally, for ge^ — 0, DM 
constraints apply essentially to the product g^g^j^. The 
individual values matter only if im^ ^ Mu, in which 
case the value of the decay width Tjj is relevant. 

Only a narrow range of couplings is allowed in Fig. 1 
on account of an S*— wave contribution surviving for our 
choice of gei, = 0. Even though proportional to rrig, for 
small m-^ it would lead to a present annihilation cross sec- 
tion above the range Q, unless 2m^ is quite close to Mu. 
In the latter case, the kinetic energy of the DM particles 
in the early universe was sufficient to allow efficient an- 
nihilation through the exchange of on-shell [/—bosons, 
whereas those at the galactic center are so slow that 
they can only annihilate through the exchange of off- 
shell U~s. This also gives the required large enhance- 
ment of the cross section at decoupling relative to that 
at the galactic center, even for a non-vanishing S*— wave 
contribution. However, if 2171^ is too close to (but still 
below) Mu, the relic density constraint will require 
very small couplings, too small to satisfy the constraint 
(O. For small Mu, the allowed range of values is, 
therefore, very narrow, e.g. 0.72 MeV < < 0.76 MeV 
for AIu ~ 2 MeV. For larger Mu, and correspondingly 
larger m^, the S~ wave contribution oc ml becomes less 
important, and a wider range of values of is allowed, 
e.g. > 11 MeV for Mu = 60 MeV; this lower bound 
on increases only slowly for even larger values of Mu ■ 
Recall that values of close to the upper bound of 20 
MeV are preferred 01 since they allow to describe an ex- 
cess of MeV photons from the galactic center. 

With the S*— wave contribution vanishing for a scalar x, 
the allowed parameter space is much wider (Fig. 2). The 
DM constraints are now compatible with the entire range 
of rriy, except for values very close to Mu /2 where today's 
Cann comes out too Small (for 2 my just below Mu) or too 
large (for 2m^ just above Mu) |2l| if the couplings are 
chosen to satisfy the constraint For example, for 

Mu = 2 MeV, the range 0.91 MeV < < 1.04 MeV is 
excluded. Note that in both Figures, smaller couplings 
correspond to larger (smaller) if 2m^ < {>)Mu. 

Since it is natural to assume g,yj^ = for our choice of 
ge^ = 0, the only relevant model-independent laboratory 
constraints come from processes involving only electrons, 
the most sensitive being the anomalous magnetic mo- 
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FIG. 2: Parameter space of the model with a complex scalar 
as MeV Dark Matter x annihilating through the exchange of 
spin-1 U bosons, for g^-^ = = and g^ = 1- Notation is as 
in Fig. 1, except that the indicated sensitivities are now those 
that can be achieved at the B— factories. 



ment of the electron Using the analytical results of 
[l^ and comparing with the most recent results for SM 
prediction and measurement we find 
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at 95% C.L. The resulting upper limit on the prod- 
uct f/e^ is shown as dotted (blue) line in Figs. 1 
and 2. While the data weakly favor a small positive 
contribution^^, [/—boson loops can account for it only 
if 5ei, — ffejj • This would either imply a sizable U vv cou- 
pling, or a model where the U gauge group is embedded 
non-trivially into the electroweak gauge group, render- 
ing the construction of a complete, renormalizable model 
more complicated. As long as ^ei, is zero, {ge — 2) imposes 
a rather severe constraint (Figs.l&2) entirely indepen- 
dent of g^^^. With DM-constraints operating essentially on 
the product gx9eR i the parameter space that also satisfies 
the {ge — 2) constraint becomes larger for larger g^. To- 
gether, the constraints exclude scenarios with g^ <^ gen- 
In the opposite limit, namely g^ = 1 (3> gen), they to- 
gether exclude scenarios with Mjj much above 200 MeV, 
as shown in Fig. 2. Other constraints on the [/—boson 
have been discussed in the literature however, 
they are more model-dependent in that they involve cou- 
plings that are not required from DM phenomenology. 
Hence we ignore these constraints here, and instead turn 
to a discussion how this model can be tested at e+e^ 
colliders. Scenarios with g^ = but ge 7^ may also be 
favored by BBN [3. 

As noted above, the U boson must couple to e+e~ pairs 
in order to explain the excess flux of 511 keV photons. 
The process e+e^ ^ U'y will therefore have a non- 
vanishing cross section if the cms energy > Mu: 

v2n 
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where a is the fine structure constant and = is the 
emerging angle of the photon. Since existing constraints 
imply Mu < 0.2 GeV, the produced U bosons can decay 
only into e"'"e~ pairs, pairs of DM particles Xj or Per- 
haps neutrinos. This gives rise to two possible collider 
signatures, with e+e~7 and 7-I- "nothing" final states 
and we explore the detectability of both. We focus on 
low-energy colliders, since the signal cross section will 
drop cx 1/s for s ^ M^. Specifically, we analyze the two 
signatures at the $ factory DA$NE, which operates at 
y/s = = 1.02 GeV, and at the B factories, which op- 
erate at y/s — 10.6 GeV; the reduced cross section at the 
latter is over-compensated by the higher accumulated lu- 
minosity (~ 500 fb^^ for both S— factories combined, as 
compared to ~ 500 pb^^ at the $ factory). 

The e"'"e~7 final state receives a large contribution 
from 0{a^) QED processes. But whereas the signal 
events have invariant mass of the outgoing e+e~ pair 
very close to Mjj, the background distribution has 
peaks at Mee — \fs (from i— channel diagrams with soft 
7 emission) and at a few me (from s— channel diagrams). 
We therefore require that (i) the produced particles must 
not be too close to the beam pipe, | cos 6*4! < 0.9 for 
i = e±, 7; (ii) Mee S [M^-l MeV, Mu + l MeV], where 
the spread is given by the mass resolution of the KLOE 
detector |19]; we assumed that the BaBar and BELLE 
detectors at the B factories have similar resolution. The 
second cut implies that the photon is quite energetic, 
since — {s — M^^)/{2y/s). The signal is considered 
detectable if A'^signai > 5^iVbckgd- 

The resulting sensitivity limits are indicated by the 
(upper) dashed (dark green) curves in Figs. 1 and 2. 
Signal and background have been calculated with the 
CompHEP package augmented to include U bosons. 
Note that these are maximal sensitivities in that we as- 



sume a branching ratio B{U 



) = 1. This may 



well be realistic for 2my^ > Af[/, but for the assumptions 
made in the Figures is not realistic otherwise; since we as- 
sumed g^ ^ gejj, , ffei, in these plots, the invisible U XX 
decay mode will dominate if it is open. However, since 
the e~'"e~7 final state is background dominated, the sen- 



sitivity limit on ^e^ only scales like B(U 



0- 



-0.25 



(6) 



Thus, DA<i>NE can probe couplings gen down to about 
10"'^, whereas the B factories would be sensitive to cou- 
plings as small as 3 • 10^''. In both cases the sensitiv- 
ity gets somewhat worse at small Mjj as the background 
peaks at small Mee. Note that these sensitivity limits are 
completely independent of g^, and of the nature of the 
DM particle (scalar or Majorana fermion). In particular, 
for the scenario of Fig. 1, the B factories should be able to 
probe the entire parameter space with 2m^ > Mjj in this 
channel. The scenario in Fig. 2 is the worst case scenario 
for collider experiments. With the DM constraints essen- 
tially fixing the product of couplings geadx (for .9ei, — 0), 
a large g^, as in Fig. 2, therefore leads to small ge^, and 
hence small production cross sections. 
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As already noted, for g-,^ ^ we expect a large invis- 
ible branching ratio for the U boson if 2m5^ < Mjj. The 
signal, then, consists of a single monochromatic photon 
with Ej = {s — M^)/(2a/s). Unfortunately, the experi- 
mental resolution on is considerably worse than that 
for Mee 0. On the other hand, the physics (SM) back- 
ground now comes from vD^ final states, and is thus 
0{aGps). After simple acceptance cuts, E^ > 100 
MeV, |cos0| < 0.9, the background is already com- 
pletely negligible at the <i> factory. Even at the B fac- 
tories we expect <C 1 background event once we require 
|cos6'| < 0.9, E-^ > 0.5Vi- 200 MeV ■22]. In other 
words, the photon plus "nothing" signal is rate, rather 
than background, limited. We neglect instrumental back- 
grounds here as these are very specific to the experiment. 

The corresponding sensitivity limits are shown by the 
(lower) dashed (maroon) lines in Figs. 1 and 2. We 
again show the maximal sensitivity, i.e. here we assumed 
100% branching ratio for invisible U decays. We see 
that DA$NE can probe a coupling > 8 • 10~^ in 
this channel, while the B factories would be sensitive to 
gen ~ 3 • 10~^, independent of g^. In particular, the 
B factories would probe the entire parameter space of 
Fig. 1 with Mij > 4 MeV in this channel. Even in the 
worst case scenario of Fig. 2, much of the DM-allowed 
parameter space would lead to a detectable signal. 

A dominant invisible decay mode for the U is not real- 
istic if 2m^ > Mjj , unless U also has significant couplings 
to neutrinos, which however is disfavored. We therefore 
also investigated XXI production through off-shell U ex- 
change. The signal cross section is now proportional to 
the product genQx (f^i' 9eL — 0), just like the DM annihi- 
lation cross section © . Using a modified CompHEP, we 
find a detectable signal at the i?-factories (> 10 events 
with \ cose^\ < 0,9, E-^ > 0.5^5 -200 MeV in 500 fb-i) 
if gen.gx ^ 2.4 • 10^^. This limit depends only weakly on 
Mjj and m^, so long as rriy^ is not very close to Mjj/2. 
This would be sufficient to probe at least the upper end 
of the DM allowed region with 2to^ > Mjj in Fig. 2. 

In summary, we carefully delineated the allowed pa- 
rameter space of models with MeV Dark Matter whose 
annihilation is mediated by the exchange of a spin-1 U 
boson. Model parameters must be chosen such that the 
correct thermal relic density and the correct present an- 
nihilation cross section are reproduced; the latter is moti- 
vated by the signal of 511 keV photons from the center of 
our galaxy. The parameter space is further constrained 
by the anomalous magnetic moment of the electron. We 
found that these models can be tested decisively at exist- 
ing low-energy e~^e~ colliders if the U boson has similar 
(or greater) coupling strength to electrons as to DM par- 
ticles (Fig. 1). Such models would be relatively easy to 
construct by introducing an additional gauge group with 
small coupling constant. Models where the coupling of 
the U boson to DM particle is much stronger than that 



to electrons are much more difficult to probe at collid- 
ers (Fig. 2); such a pattern could emerge if U couples to 
electrons only through mixing with SM gauge bosons, but 
has direct coupling to DM particles. The single photon 
plus "nothing" channel allows to probe much of the pa- 
rameter space compatible with the DM constraints even 
in this unfavorable situation. We therefore come to the 
rather surprising conclusion that the solution of the Dark 
Matter puzzle might be found at the existing $ and B 
factories. 
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